Geometric quantum computation using fictitious spin-1/2 subspaces of strongly 

dipolar coupled nuclear spins 
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Geometric pliases liave been used in NMR, to implement controlled phase shift gates for quantum 
information processing, only in weakly coupled systems in which the individual spins can be identified 
as qubits. In this work, we implement controlled phase shift gates in strongly coupled systems, by 
using non-adiabatic geometric phases, obtained by evolving the magnetization of fictitious spin-1/2 
subspaces, over a closed loop on the Bloch sphere. The dynamical phase accumulated during the 
evolution of the subspaces, is refocused by a spin echo pulse sequence and by setting the delay of 
transition selective pulses such that the evolution under the homonuclear coupling makes a complete 
2n rotation. A detailed theoretical explanation of non-adiabatic geometric phases in NMR is given, 
by using single transition operators. Controlled phase shift gates, two qubit Deutsch-Jozsa algorithm 
and parity algorithm in a qubit-qutrit system have been implemented in various strongly dipolar 
coupled systems obtained by orienting the molecules in liquid crystal media. 

I. INTRODUCTION 



The concept of using quantum systems for information processing was first introduced by BeniofF Qj. In 1985 
• I— I 

Deutsch described quantum computers which exploit the superposition of multi particle states, thereby achieving 
massive parallelism ^j. Researchers have also studied the possibility of solving certain types of problems more 
efficiently than can be done on conventional computers y, |j, |5{ . Quantum computing and quantum information 
processing, requires the ability to execute conditional dynamics between two qubits P| . Several techniques are being 
exploited including nuclear magnetic resonance The experimental realization has been limited by 

many factors, such as decoherence . d ifficulty in achieving higher number of qubits and controlled manipulation of 
qubits with higher fidelity [ij, |l3|, |lj, |l5|, |la, llTl. A recent development in this field is to use geometric phases to 
build fault tolerant controlled phase shift gates [ig • 
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Berry's discovery of geometric phase Q| accompanying cyclic adiabatic evolution has triggered an immense effects in 
holonomy effects in quantum mechanics and has led to many generalizations. Simon explained that this geometric 
phase could be viewed as a consequence of parallel transport in a curved space appropriate to the quantum system. 
In nuclear magnetic resonance, geometric phase was first verified by Pines et.al. in adiabatic regime |2l| . A similar 
approach was used by Jones, to implement controlled phase shift gates by geometric phases in a two qubit system 
formed by a weak J couplin g [221 123| . However the adiabatic condition is not satisfied in many realistic cases because 
of the long operation time 24] . Hence it is difficult to experimentally realize quantum computation with adiabatic 
evolutions, particularly for systems having short decoherence time. To overcome this disadvantage, it was proposed 
to use the non- adiabatic cyclic geometric phase (Aharonov and Anandan phase 26]) to construct geometric quantum 
gates [2^. These gates not only have faster gate operation time, but also have intrinsic features of the geometric 
phase, and hence robust against certain types of operational errors For a non-adiabatic cyclic evolution, 

the total phase difference between the initial and final states, consists of both the geometric and dynamical phases. 
Therefore, to obtain only the non-adiabatic geometric phase, one has to remove the dynamical component. Non- 
adiabatic geometric phase in NMR was also first verified by Pines et.al Non adiabatic geometric phases in NMR, 
were used to implement controlled phase shift gates, Deutsch-Jozsa (DJ) and Grover algorithms in a two qubit system 
formed by weak J coupling [3| • 

Most of the NMR Quantum information pro cessing (QIP) experiments have utilized systems having indirect spin- 
spin couplings (scalar J couplings) |7l Isl M. llOl l29l|. Since these couplings are mediated via the covalent bonds, the 
number of coupled spins and hence the number of qubits is limited to a few aubits Another approach is to use 
direct dipole-dipole couplings between the spins which are larger in magnitude [sil Isli Is^l Is^ Is^ Is^ . In hquids the 
dipolar couplings are averaged to zero due to rapid isotropic reorientations of the molecules and in rigid solids there are 
too many couplings yielding broad lines 39]. In molecules partially oriented in anisotropic media, like liquid crystals, 
one obtains partially averaged intra molecular dipolar couplings with only a finite number of dipolar coupled spins, 
yielding finite number of sharp NMR resonances 30] . However in such cases often the dipolar couplings are large or 
comparable to chemical shifts differences between the coupled spins, leading to spins becoming strongly coupled. The 
method of spin selective pulses and J evolution used in NMR-QIP of weakly coupled spins (liquid state NMR-QIP), 
leads to complications in strongly coupled spins because of the non secular terms in the Hamiltonian. Therefore such 
systems have been utilized for quantum information processing, by using transition selective pulses and treating the 
2" non-degenerate energy levels as an "n" qubit system. Preparation of pseudopure states, one qubit DJ algorithm 
and quantum gates which do not need control of phases, have been implemented in such systems using transition 
selective pulses [s^il^ls^. Recently an 8 and a 12 qubit systems have been achieved using dipolar coupled spins in 
molecules oriented in liquid crystals 0, [s^ . 
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In this paper, we implement controlled phase shift gates, two qubit DJ algorithm and the parity algorithm in 
strongly dipolar coupled spins, by using non-adiabatic geometric phases, obtained by evolving the magnetization of 
fictitious spin-1/2 subspaces (formed by pairs of energy levels). This method requires efficient transition selective 
pulses and refocusing of Hamiltonian evolution. The physical systems chosen, contain a single non secular term h-Ij 
in the homonuclear dipolar Hamiltonian Dij(31zilzj — h-Ij), which is refocused by setting the delay of a transition 
selective pulse such that the evolution due to the coupling makes a complete 2tt rotation. To the best of our knowledge, 
this is the first implementation of quantum algorithms in strongly coupled nuclear spins by using geometric phases. 

In section(II), we outline the non-adiabatic Geometric phases by using both, Bloch sphere approach and single 
transition operators In section(III), we describe the implementation of controlled phase shift gates in ^^CH^I 
partially oriented in ZLI-1132 liquid crystal. In section(IV) two qubit DJ algorithm is implemented in dipolar coupled 
equivalent protons of oriented CHj,CN . In section (V) parity algorithm is implemented in a qubit-qutrit system of 
CH2FCN partially oriented in ZLI-1132, and the results are concluded in section(VI). 



II. NON-ADIABATIC GEOMETRIC PHASES IN NMR 



In order to describe the non-adiabatic geometric phases in NMR, we discuss the evolution of the magnetization of 
sing le quantum transitions, under transition selective pulses, using the approach of single transition operator algebra 
|40j, briefly outlined below. 

Consider a 4 level system (fig. la), consisting of 6 pairs (r,s) of energy levels (r=l,2,3, and r < s < 4 ). Each pair 
of energy levels can be considered as belonging to a fictitious spin-1/2 subspace (virtual two level system), and the 
dynamics of the magnetization of the system can be explained by using single transition operators |4c| . Here we 
assume that the single quantum magnetization of the subspaces (1,2), (2,3) and (3,4), is directly observable, hence 
one can easily manipulate those subspaces. The angular momentum operators (single transition Cartesian operators) 
la'''^ (Q!=x,y,z) have been defined as 

W'^'%)^\{5^r5,r-5^s5,s), (1) 

where i, j = 1, 2, 3 and 4. From this definition it follows that J4fl|, 

= I^^^'-), li^-'^) = -/^^■'■) and I'f''^ = -4''"^- (2) 



In the following, geometric phases are explained by using Bloch sphere approach Later, we describe the 

commutation relations of single transition operators eqn.J^I, and show that geometric phases can be obtained by 
using pair of transition selective tt pulses, with suitable phases in the xy plane. The removal of dynamical phase 
(phase acquired due to internal Hamiltonian evolution) is discussed in sections III, IV and V. 

From equation JQ), one can see that, the angular momentum operators of (r,s) subspace, after removing rows and 
columns containing zeroes, are Pauli spin matrices and hence the basis states of (r,s) subspace can be considered as 
|0) and which can be represented on a Bloch sphere fig.^D). In this paper we follow the convention that the 
lower state 'r', in the (r,s) subspace, represents |0) and the upper state 's' represents In general, every two level 
sub-system at equilibrium, can be considered to be in a pseudopure state QH. Hence the (r,s) subsystem can be 
considered as in the pseudopure state |0). In such systems, to encode the phase information, we create initial states, 
consisting of superposition states, in each of the subspaces (1,2), (2,3) and (3,4), given by. 



IV'rs) = — ;=(|0) + |l))rs, where r = l,2,3 and s = r + 1. 
v2 



(3) 



The state of (r,s) subspace, IV'rs) (eqn|2Jl, consist of coherences llf''^'^ (r=l,2 and 3, s=r-|-l ), which can be created 
experimentally by applying a (j: /2)y pulse on the equilibrium state of the whole system. The coherence /i'^'*^ can in 
general be written as. 



(4) 



where l'^^''^^ is the identity matrix of (r,s) subspace and \'iprs){'4'rs\ represents the density matrix cr^r.s)- Since the 
identity matrix does not contribute to NMR signal, the density matrix (J(r,s) = |'0rs)('0rs| represents I^'"^^ and vice 



versa. 



The state ji/irs), can be manipulated, such that it goes through a closed loop on the Bloch sphere, thereby acquiring 
a geometric phase (fign^). The geometric phase acquired is respectively e*'^/^ and e^'^/^ for the states |0) and 
where 51 is the solid angle subtended by the closed loop at the center of the Bloch sphere (figHs) [2^ 12^- The unitary 
operator for the above operation is a phase shift gate, given for example, for the (2,3) subspace as. 



C/(2,3) 



/I 

e*"/2 



Vo 



0\ 


1/ 



(5) 
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The state \iprs), under the above unitary operator U(2,3), is transformed to \iprs), given by, 

|V'^3) = i=(e'^V2|o) + e-fV2|i))^3 

|^^4) = i=(e-"/2|0) + |l))34. (6) 
The states |?/'rs)i iu eqnlHI can be written in the density matrix form [(t^^ = |V'rs)(^rsl]j as, 

1 / 1 e-'^y^\ I f I e'^\ 1 / 1 e-'"/2\ 

and 0-^4 = o ,o /o _ ■ (7) 
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In the foUowing we show that, the geometric phase gate U(2,3) can be obtained by, cychc evolution of (state 
|V'23))7 by means of selective tt pulses on transition (2,3). 

The effect of transition selective pulses can be explained by using the commutators of single transition operators 
given in eqn.(Q. If the coherence and the r.f. pulse involves the same transition, then the transformation follows the 
usual commutation relation [40I 1 , 

[I^^''\I^;''h=^4'^'^\ (8) 

where (a,/3,7) is a cyclic permutation of (x,y,z). 

The commutator in eqn|Sl implies the property that 

which can be interpreted as, rotation of /i^'"'' towards around the 7 axis of (r,s) subspace. 

If the pulse is applied on the transition (r,s), and a connected transition (t,r) is observed, then the commutation 
relations become j40j, 

l-^X T^X J 2 y 

and [4*''-),4*-)] = ^/^^^ (10) 
The commutator relations of eqn.(lO) imply that: 

and e-^''^^°'/(*^'^)e^^^^"' = 4*'''^cos(^) + I^'-'\sin{^-). (11) 
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The first commutator in eqn.l^UJ, implies that, the coherence Ix*'^^ is rotated towards /i*'*'', when the pulse is 
applied on transition (r,s), about y-axis. Whereas the second commutator in ean. (|10|l . implies that, ly^'^^ is rotated 
towards Ix''^\ under the r.f pulse, applied on transition (r,s) about the x-axis. Unlike in eqn.®, the factor (1/2) 
in ean. Hll)|) . indicates that, the angle of rotation is halved, when the transformed operator and the rotation operator 
have only one index in common. Furthermore the commutation relations of ean. l|lU|) are not cyclic. The commutation 
rules, given in ean. HlU|) . can also be visualized in three dimensional subspaces spanned by three orthogonal operators 
each, as shown in fig|2 for the subspaces (1,2) and (3,4), when the pulse is applied on (2,3) subspace. It may be noted 
that, in certain rotations the single quantum magnetization will evolve into a double or zero quantum coherence. 

We now apply two consecutive selective tt pulses on coherence Ix ' (or state \1p23), eqn© and &g-^p)), with 
phases respectively given by, y and {y + ir + cj)) in the xy plane. The evolution of the coherences /i^'^^ and /^'^'^^ 

can be calculated by using the commutation relations given in eqn.s ^ and H1U|I (figlU, 
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cos(20)/F'^^ - sin{2(j))4^'^\ 
cos4>l'i''^^ + sin(f)Iy^'^\ 

The unitary operator associated with above operation (eqn.l2) is identical to the phase shift gate U(2,3) of eqn|31 
with fl/2 ~ <j). In other words, the transformed coherences of equation (12), are identical to {a'^^ — (1/2)/'^''''')) 
of eqnEI with = fl/2. Hence one can conclude that, two selective tt pulses, having phases y and (y + tt + 0) 
respectively, applied on coherence , evolves the state |'023)7 over a closed loop on the Bloch sphere, and the solid 
angle subtended by the closed loop, at the center of the sphere is f2 = 2(j>. However, in case of a 7r-phase shifted gate 
{(j> — tt), the two TT pulses, (7r)y and (7r)(y^^^7r=i/)7 shown in fig.lQb), can be combined in to a single (27r)j, pulse. Such 
pulses are used in sections HI, IV and V. 

It is to be noted that, if the coherence of the (2,3) subspace is in the xy plane, then the phases of tt pulses 

will be {y + 6') and {y + 9' + tt + </)), and a 7r-phase shifted gate is achieved by applying (27r)(j,^g/) pulse. And if the 

(2, 3) 

(2,3) subspace is at equilibrium state, Iz , the phases of tt pulses will be {9") and {9" + tt + (p), where 9" can be 
any axis in the xy plane, and a 7r-phase shifted gate is achieved by applying (27r)e" pulse. 

From the above discussion it is concluded that, cyclic evolution of a two level subspace, by means of transition 
selective tt pulses, acquires a geometric phase which also effects the other connected subspaces. 
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III. EXPERIMENTAL IMPLEMENTATION OF THE CONTROLLED PHASE SHIFT GATE 

The system chosen is, dipolar coupled carbon and protons of the methyl group of ^^CH^I, partially oriented in 
ZLI-1132 liquid crystal. The three protons are, chemically and magnetically equivalent, and the total Hamiltonian is 
given by. 



ij = l,2,3, i<j, 

where coh and uc are chemical shifts of protons and carbon respectively, Dhh is the residual dipolar coupling 
between identical protons (homonuclear coupling), and Jch and Dch are couplings of carbon (I) with protons (S). 
The scalar J coupling between identical protons does not effect the spectra in this case |39| . 

The energy level diagram (figO)) of this system consists of 16 energy levels, eight of which belong to a symmetric 
manifold, and the other eight to two asymmetric manifolds The eight energy levels of the symmetric manifold 
can be considered as a 3-qubit system. In this system recently pair of pseudo pure states (POPS) and quantum 
gates (other than controlled phase shift gates) have been implemented by applying transition selective pulses [s^l^. 
In this work we implement controlled phase shift gates by using non-adiabatic geometric phases obtained by using 
transition selective pulses. 

The equilibrium and ^^C spectra are given in fig.l@Ji) and fig.Qa) respectively, wherein the various transitions 
are labeled in accordance with the notation used in fig|2l Since several transitions of symmetric and asymmetric 
manifolds, are degenerate, selection of symmetric states (SOSS) is performed to select the magnetization of the 
symmetric manifold |3l| . However, since in the present experiments only the proton transitions are observed, we 
perform only a partial SOSS, in which the proton magnetization of symmetric manifold, is selectively retained, as 
explained below. 

The entire experiment consists of 4 steps (as shown in fig|3J). 

(i) Selection of symmetric states (SOSS): The proton magnetization from asymmetric manifolds is saturated by 
using the pulse sequence {7T)^''-''^'^-Gz-{'n')'^^''^* (figlp , where (7r)''3.''4 represents two selective tt pulses on and 
transiti ons, and Gz is the pulsed field gradient |3l| . The tt pulse on transitions and acts as a 7r/2 pulse on h'^ 
and /14 |32l |. and a gradient dephases the coherences /13 and h'^. The second tt pulse brings back the magnetization 
of symmetric manifold to its equilibrium value. Proton spectrum obtained by a small angle pulse followed by SOSS 
(figEfc)) has the intensities in the expected ratio 3:3:4:4:3:3, confirms the retention of proton magnetization only from 
the symmetric manifold. 

(ii) An on resonance {■K/2)y pulse is applied on the proton channel, which creates the proton coherences, hi, /i2, 
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hs, /is and /ig . These coherences does not evolve under the chemical shift, since pulse is applied at resonance. 

(iii) Controlled phase shift gates: As described in the previous section, cyclic evolution of a two level subspace by 
means of transition selective tt pulses, acquire a geometric phase. In order to get only the geometric phase, one has 
to refocus the dynamical phase acquired due to internal Hamiltonian evolution. Spin echo sequence (t — vr — r) is 
applied on protons (figEJ to refocus hetero nuclear {^^C-^H) coupling (dipolar and J). The geometric phase gates are 
obtained by applying two selective tt pulses on carbon transitions, (figE)), and the duration of the pulses is chosen 
such that the evolution of proton coherences under the homonuclcar coupling makes a complete rotation (in multiple 
of 27r). 

(iv) Detection: The phases of proton coherences hi...he (figures El and [7|), are only due to cyclic geometric evolution 
of subspaces, confirming the implementation of phase shift gates, as explained below. 

In fig.®, geometric phases are obtained in (1,2) subspace , by applying two consecutive phase shifted selective tt 
pulses on transition Ci. The phase of the first pulse is y, and that of second pulse is {y + n + 0), where 4>—7t/2, tt, 
37r/2 and 27r, for fig.sE^a), Efb), E^c), andEl^d) respectively. The unitary operator obtained by the above operation 
is, U(l,2)=diag[e"^,e~*'^,l,l,l,l,l,l ]. As expected, only the proton coherences hi and h2 are effected by the unitary 
operator U(l,2), whereas the other proton coherences remain un affected. For example in figllHt)? the coherences 
Ix'''^\hi) and li'^''^\h2), are transformed to ly^'^'' and — /^^'^^ giving dispersive signals of opposite sign, confirming 
the unitary operator U(1.2), for (j> = 7t/2. 

Figure(|7I) shows the implementation of (7r)-phase shift gates on two unconnected subspaces; U[(3,4),(5,6)]=diag[l,l,- 
1,-1,-1,-1,1,1 ], U[[(l,2),(3,4)]]=diag[-l,-l,-l,-l,l,l,l,l] and U[(l,2),(5,6)]=diag[-1,-1, 1,1,-1,-1, 1,1 ] (fig.sdi, andCt 
respectively). Each of these gates require two (27r)j, pulses on unconnected subspaces. U[(3,4),(5,6)] is implemented by 
applying {2Tr)y pulses on transitions C2 and C3, U[(l,2),(3,4)] is implemented by applying (27r)j, pulses on transitions 
Ci and C2. Whereas U[(l,2),(5,6)] is implemented by applying (27r)j^ pulses on transitions Ci and C3. The signs 
of various transitions in the observed spectra of fig. El confirm the implementation of these phase shift gates. For 
example, the U[(3,4),(5,6)] gate inverts only the coherences hi, h2, h^ and h^. 



IV. IMPLEMENTATION OF DEUTSCH-JOZSA (DJ) ALGORITHM 

D J algorithm determines in a single query, whether a given function is constant or balanced [isi . A function 
is "constant" if it gives the same output for all inputs, and "balanced" if it gives one output for half the number of 
inputs and another for the remaining half. Classically for an n bit binary function, at least (2"^^ -I- 1) queries are 
needed to determine whether the function is constant or balanced, whereas the DJ algorithm requires only a single 
query. The Cleve version 44] of the DJ algorithm requires an extra qubit (ancilla qubit) and uses controUed-not 
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gates, whereas Collins version does not require the ancilla qubit but needs controlled phase shift gates l45l . While 
the Cleve version of DJ algorithm has been implemented both in weakly and strongly coupled systems 

013: the 

Collins version has been implemented only in weakly coupled systems l43!|. Here we implement the Collins version 
of the DJ algorithm on a strongly coupled two qubit system, formed by dipolar coupled methyl protons of CH3CN, 
oriented in ZLI-1132 liquid crystal. For a two qubit system there are 2 constant and 6 balanced fmictions, table(l). 
The quantum circuit of the Collins version of the DJ algorithm is shown in figlHl The algorithm (figlHl) starts with 
a pure state (pseudo pure state in NMR) |00) which is converted to a equal superposition state of the basis states 
|00), |01), |10), and by applying a pseudo Hadamard gate on both the qubits. Thereafter a unitary operator f7/. 
(controlled phase shift gate) is applied, followed by detection. Theoretically, after Uf- one has to apply the Hadamard 
gate. In NMR, a Hadamard gate is replaced by a pseudo Hadamard gate, which is implemented by a (7r/2)j, pulse, 
and for detection another (7r/2)_j, pulse is needed. These two pulses cancel each other and the result of the algorithm 
is available immediately after Uf-. The Unitary operator Uf- can be written as, 



Uf. 



^(_l)/>(oo) 

(-l)/.(oi) 

(-i)/.(io) 

V 



(14) 



(-!)/'(") / 

where i=l and 2 represent the "constant" function and 1=3,4,5,6,7 and 8 represents the "balanced" functions, table 



The unitary operators, J7/j, Uf^, Uf^ and Uf.^ are identical, up to an overall phase factor, of e*'^, to Uf^, Uf^, Uf^ 
and Ufg respectively. Hence we implement only the unitary operators Uf-^ , Uf^ , Uf^ and J/y.. . 

The schematic energy level diagram of the dipolar coupled methyl protons of CH^CN (without C-13 labeling) 
oriented in ZLI-1132 liquid crystal, is given in (figlHll- Transitions 1, 2 and 3 of symmetric manifold gives rise to 
three single quantum transitions at frequencies lOq — 3-D, lOq and uJq + 3Z?, with the intensity ratio 3:4:3, where D is 
the partially averaged (residual) dipolar coupling between the protons of CH^CN and lUo is their Larmor frequency. 
Transitions of asymmetric manifold also have the frequency lUo identical to transition 2. The resultant spectrum is 
a 1:2:1 triplet as shown in fia llUh . The symmetric manifold of this system, consisting of four energy levels, can be 
treated as a two qubit system, and the labeling is given in figlHl 

The total experiment consists of 4 steps (as shown in fig. Ill|l . 

(i)Selection of symmetric manifold: The magnetization of asymmetric manifolds which contribute to central tran- 
sition, is saturated by applying a n pulse on transitions 2, which acts as a n/2 pulse on transitions 4 and 5, a 
gradient pulse dephases the coherences, and the magnetization of the symmetric manifold is retained to its equilib- 
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rium value, with another tt pulse. The relative integrated intensities of figllUb (in the ratio 3:4:3) confirm the selection 
of symmetric manifold. 

(ii) Creation of pseudo pure state (PPS) |00) [s^: The pseudo pure state |00) is prepared by applying a tt pulse on 
transition 3 (fig. Ill() , followed by a 7r/2 pulse on transition 2, a gradient pulse is applied to daphase the coherences. 
The proton spectrum of figllUb. confirm the creation of the |00) PPS. 

(iii) Superposition state followed by Uf^: After the creation of |00) PPS, a pseudo Hadamard gate ( {n/2)y on 
resonance pulse) creates a superposition of 4 basis states, and due to on-resonance pulse the coherences do not evolve 
under the chemical shift. It may be noted that unlike in weakly coupled spins, the state created here is not in uniform 
superposition of eigenstates, since the coefficients of various eigenstates are different IsiJ. However, as is shown here, 
the created coherent superposition can be utilized for quantum parallelism, to distinguish different classes of functions. 
The state of the system after a (7r/2)y pulse is [V') = (l/2y2)(|00) + %/3|01) + V3|10) + |11)). 

The unitary operator Uf^ of ean. (|14(l . is implemented by using non-adiabatic geometric phases, obtained by tran- 
sition selective {2-K)y pulses. The duration of transition selective pulses (fig. 111(1 is set such that the evolution due to 
homonuclear — ^ H dipolar coupling, makes a complete 27r rotation. U f-^ is a unit matrix which require no pulses. 
U and are implemented by applying a single transition selective (27r)j^ pulse on transitions 2 and 3 respectively. 
Uf^ is implemented by applying two transition selective (27r)y pulses on transitions 2 and 3. The final state 
(=[//. IV)) can be written as. 
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(15) 

(iv) Detection: The phases of coherences 1, 2 and 3 confirm the final states IV/J- For constant function (IV/i)), 
none of the peaks are inverted ffigll2b'). whereas for the balanced fimctions (IV/a), iV/s) a-nd iV/?)), atleast one of 
the peaks is inverted ffigll2btT^ll2ti respectively). The spectra of figure El confirm the implementation of Collins 
version of two-qbit DJ algorithm in the dipolar coupled spins by using non-adiabatic geometric phase gates. 



V. IMPLEMENTATION OF PARITY ALGORITHM 



The parity of a binary string X={xi, , Iat}, where Xi is or 1, is even (odd) if x consists of even number 

(odd number) of I's. The parity of X can be written as a Boolean function P(X)=a;i® ©xjv . A classical computer 
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has to call the oracle with each of the N possible inputs to determine P(X). Beals et. al and Farhi et. al showed 
that in a quantum computer the minimum number of oracle calls is N/2 |47l l48l| . The parity algorithm has been 
implemented by Suter et.al on a weakly coupled two qubit system. They have also shown that in case of ensemble 



quantum computers the complexity is further reduced |49j . In the quantum version, the parity of the binary string X 

is encoded in the Black box (oracle), such that the black box consists of N boolean variables, X~{xi, 

The task is to determine the parity of the binary string, that is encoded in the oracle. 

Recently it has been demonstrated that qutrits can be useful for certain purposes of quantum simulation, quantum 
computations and quantum communication, a qutrit has three states |0), |1) and |2) |5Cl ISlL l52l l53l 154. ISSj . In this 
work, we use a quantum circuit similar to the one used in ref. 49] and implement a qubit-qutrit parity algorithm in a 
single iteration of the oracle, starting with a mixed state. Non-adiabatic geometric phase gates are used to implement 
the required oracles. The qubit-qutrit system used here is obtained by orienting flouro-acetonitrile in ZLI-1132 liquid 
crystal. The two identical protons have one asymmetric and three symmetric eigen states (labeled as 0, 1, 2 in figll4|l. 
The symmetric manifold of protons can be treated as a qutrit and the flourine nucleus as a qubit, yielding a 6 energy 
level system shown in fig. (|14|l . 

In the quantum circuit of the algorithm, the first register ffig ll3|l is assigned to flourine (qubit) and the second 
to protons (qutrit). The algorithm starts with a mixed state (ensemble of 6 basis states of a qubit-qutrit system), 
a pseudo Hadamard gate on the qubit register, creates single quantum coherences, whose phases are modulated by 
the oracle. Detection is done on the qubit register, which measures the single quantum coherences. After the pseudo 
Hadamard gate, the polarization of the qubit register consists the sum of three single quantum coherences, which can 
be written as. 



ir"' = |0)(1||0)(0| + |0)(1||1)(1| + |0)(1||2)(2|, (16) 

where the first part of each term indicates the single quantum coherence of the qubit and second part indicates the 
state of qutrit. 

The oracle in ffigll3|l is defined as 
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Table(2) shows the unitary operators of the oracle corresponding to different strings, X={xi, xg}, of various 
parities. Here, one should note that the oracles, for parities 0, 1 and 2 are identical to that of parities 6, 5 and 4 
respectively, up to an overall phase factor e'''^{— —1). Hence we implement only the oracles of parities 0, 1, 2 and 
3. For even parity, even number of terms of ean. (|16|l acquire a phase factor e*'^, and for odd parity odd number of 
terms acquire a phase factor e*'^. Hence, by measuring single quantum coherences, one can determine the parity that 
is encoded to the 6 digit binary string. 

The equilibrium spectra of and ^^F are shown in fig llSb , and 115b respectively. In the flourine spectrum, the 
intensity of central transition is twice that of satellite transitions, due to the overlap of the transition of the asymmetric 
manifold. 

The experiment consists of 3 steps (as shown in fig ll6|l 

(i) Selection of symmetric manifold (SOSS): The magnetization from asymmetric manifold is saturated by using 
the pulse sequence ffig ll6|l . known as SOSS. A gradient pulse following a tt/2 pulse on flourine, saturates the flourine 
magnetization. The flourine magnetization of symmetric manifold is retrieved by partial polarization transfer from 

to ^^F by applying selective tt pulses on H2 and transitions followed by a selective Stt/A pulse on F2 transition. 
A final gradient pulse dephases the transverse magnetization created during the process. The spectrum obtained by 
a small angle pulse after SOSS, gives the relative intensities in the expected ratio 1:1:1 ffigllSb). confirming SOSS. 
The six energy levels of the symmetric manifold can be treated as a 6 digit binary string X={xi, x^}. 

(ii) Implementation of pseudo Hadamard gate followed by oracle: A (7r/2)j, on resonance pulse, applied on ^^F, 
implements a pseudo Hadamard gate, and the created coherences does not evolve under the chemical shift, since the 
pulse is applied at resonance. The unitary operators of the oracle, are implemented by using non-adiabatic geometric 
phases. The phase shift gates d^, or (Table 3) can be implemented by applying two transition selective tt pulses 
respectively on transitions Hi or Fi with phases y and (y + tt + 0), and the duration of the pulses is set such that the 
evolution under the hetero nuclear couplings (J + 2D) makes a complete 2tt rotation. However spin echo sequence 
(r — vr — r) is also used which refocuses hetero nuclear couplings and chemical shift evolution (if any). Since the unitary 
operators of the oracle are diagonal matrices, one can sandwich different phase shift gates to get a desired phase shift 
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gate, as shown in Table(2), up to an overall phase factor. For example Of is a product of d^^^, '^lls ^ ^Jil^^ ^ '^Vi ^ '^'f^ 
(up to an overall phase factor e"/^), which is implemented by three pairs of selective tt pulses oxi^H channel and 
two pairs of tt pulses on ^^F channel, fig. 16. Similarly one can implement remaining oracles by sandwiching various 
phase shift gates, as shown in Table 2. 

(iii) Detection: Single quantum coherences of flourine are measured. Depending on the number of peaks that are 
inverted the parity is determined such that, for even parity, even number (0 or 2) of peaks are inverted, while for 
odd parity, odd number (1 or 3) of peaks are inverted. The figures [TTT a). ITTTb). llTf c) and llTr d) indicate the odd 
parity, which are obtained by implementing 0\, Of, Of and O^'"^'^^ (Table 2) respectively. Figures ITTT e). ITTf fl. 
llTt g) and llTr hl indicate the even parity, which are obtained by implementing I, 0^''^\ 0^2'^^ ^'I'i O^'^^ (Table 2) 
respectively. These spectra clearly confirm the implementation of the parity algorithm in this qubit-qutrit system, 
utilizing non-adiabatic geometric phase gates. 



VI. CONCLUSIONS 



Conventionally, controlled phase shift gates which are basic elements in many quantum circuits, are implemented 
by the combination of qubit selective pulses and evolution under the weakly coupled Hamiltonian. In order to achieve 
higher number of qubits in NMR, one explores dipolar couplings which are larger in magnitude, yielding strongly 
coupled spectra. In such systems since the Hamiltonian consists of non-secular terms it is difficult to apply qubit 
selective pulses. The present method to implement controlled phase shift gates, does not require qubit selective pulses 
and evolution under Hamiltonian. While in this paper, we have used the systems containing only one homo nuclear 
dipolar coupling, the method is being extended to systems containing more than one homo nuclear dipolar coupling, 
by using suitable refocusing schemes and by using strongly modulated pulses This method can also be applied 
to quadrapolar systems. This method thus contributes towards the use of strongly coupled spins for NMR Quantum 
computing. 
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Table 1. Constant and Balanced functions for a two qubit system 





Constant 
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X 
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Table 2. Unitary operators of the oracle, for various parities of a binary string X 



Parity 


X 


(oracle) 





{0,0,0,0,0,0} 


diag[l,l,l,l,l,l]=I 


1 
1 
1 


{0,1,0,0,0,0} 
{0,0,0,1,0,0} 
{0,0,0,0,0,1} 


Ol = diag[l, -1, 1. 1, 1. 1] = d^^'/^d^f 
Of = dtag[l, 1, 1,-1, 1, 1] = d%'.d%\d-„:^\d%'.d%' 
Ol = d^ag[l, 1, 1, 1, 1, -1] = rf^f .rf^f .d^f 


2 
2 
2 


{0,0,0,1,0,1} 
{0,1,0,0,0,1} 
{0,1,0,1,0,0} 


O^'') =rfmg[l, 1,1, -1,1,-1] = d^^ 
O^'-'^ = diag[l, -1, 1. 1, 1. -1] = d^H.-dli, 
Of'^^ ^diag[l, -1,1, -1,1,1] ^dl^ 


3 


{0,1,0,1,0,1} 


O'i^^''^ = diag[l, -1, 1, -1, 1,-1] = dli^ .d-^l'\d%\dl^ 



Table 3. Unitary operators of controlled-c?!) phase shifted gates 



d^ 

Hi{Fi) 




= diag[e^'>', e"'"^, 1, 1, 1, 1] 




= dm5[l,l,e^'^,e-''^,l,l] 




= diag[l, 1,1,1, e'"^, e"*"^] 




= dia5[e''^,l,e-''^, 1,1,1] 


4. 


= rfia5[l,e''^,l,e-''^,l,l] 


43 


= dia5[l,l,e*'^,l,e-*'^,l] 


4. 


= rfia5[l,l,l,e''^,l,e-''^] 
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FIGURE CAPTIONS 

(1) . (a) Schematic energy level diagram of a four level system, such that the single quantum magnetization of the 
subspaces (1,2), (2,3) and (3,4) is directly observable, (b) Evolution of state \ijjrs) (r,s=l,2 and 3, s=r+l), over a 
closed loop on the Bloch sphere, by applying two transition selective tt pulses on transition (r,s) with phases y and 
(y + n + n = y) . Solid angle subtended by this closed loop at the center of the sphere is, 27r. 

(2) . Rotations in subspaces spanned by orthogonal operators, of the 4 level system of fig. (la), (a), (b), (c) 
and (d) follow the commutation relations, respectively given by, li^'^"^] = -^ly^'^K [Iy"'^\ Ix"'^"^] = ■^Ix^'^\ 

ly^'^"^] = and li^''^''] — ■^ly^'^'' (eqn.lO). Each of the commutation rules follow the transformation 

property given in eqn.(ll), for example, for fig. (a) li-^'^^] — ^/y^'"^'' , means that a {(j))y pulse applied on transition 

(2,3), acts as a 0/2 pulse on transition (1,2), hence Ix is transformed to (li ' ' cos{(j)/2) + ' sm(0/2)). (Adopted 
from ref.s [39], [40]) 

(3) . The schematic energy level diagram of dipolar coupled protons and carbon of ^^CH^I, partially oriented in 
a liquid crystal. There are 16 energy levels, eight of which belong to a symmetric manifold and the other eight to 
two groups of asymmetric manifolds, hi, h[ and h'l represent proton transitions and Ci, C[ and represent carbon 
transitions. All transitions with identical subscript are degenerate, eg. C^, C'^ and C" have identical frequencies, as 
well as hi, h[ and ft,", (see fig. 4). 

(4) . Proton (a) and carbon (b) equilibrium spectra of oriented ^'^CH^I recorded at 500 MHz of proton frequency. 
The labeling of transitions are in acx'ordance with the notation used in fig. 3. The transitions /i2, /14 and /if, (and 
also hi, /i3 and /is) are split by homonuclear dipolar coupling, and the splitting is equal to '.iDuH which has a 
value 3553 Hz, fig. (4a). The C-H splitting is equal to {2Dch + Jch), which has a value 2053 Hz, fig. (4b). The 
relative experimental integrated intensities are in the ratio; for 0.98:1.00:2.00:2.02:1.00:0.99 (theoretical ratios 
are 1:1:2:2:1:1), and for ^"^C 0.98:3.02:3.00:1.00 (theoretical ratios are 1:3:3:1). (c) spectrum, obtained using a 
small angle measuring pulse, after implementation of SOSS. The SOSS is implemented here by the pulse sequence 
(^T^^h^M _ Gj, - (Tr)''^'^^ The relative experimental integrated intensities are in the ratio 2.97:3.02:4.05:4.02:3.00:2.98 
(theoretical ratios are 3:3:4:4:3:3). 

(5) . The pulse sequence for implementation of Geometric phase shift gate (GPSG) U(l,2)=diag[e"^,e~"^,l,l,l,l,l,l]. 
The pulse sequence consists of four parts; (i) SOSS:[(7r)''3''** — Gz — (tt)'*-'''*^], (n) Pseudo Hadamard gate on protons: 
{TT/2)y^ , (iii) Geometric phase shift gate: [r — tt — r], {{Tr)y^ (7r)^^^_i_0) during second t period, evolves the z 
magnetization of (1,2) subspace over a closed loop on the Bloch sphere, with solid angle 20 (iv) Detection of proton 
resonances. 

(6) . Proton spectra obtained after the implementation of Geometric phase shift gate 



18 



U(l,2)=diag[e"^,e~"^,l,l,l,l,l,l] (pulse sequence of fig. 5). (a) to (d) respectively, are for (j) = 7r/2, tt, 37r/2 
and 2n. Since the Hamiltonian evolution during transition selective pulses is rcfocuscd, the phases of the the peaks 
in the spectra are entirely due to geometric phases and confirm the implementation of geometric phase shift gates. 
The transition selective carbon pulses arc low power Gaussian pulses of duration 5.65 msec (=r/2). 

(7) . Experimental proton spectra obtained after implementation of (7r)-phase shifted gates on two unconnected 
subspaces. These gates are obtained by applying two transition selective (27r)y pulses on carbon transitions, during 
the second r period in the pulse sequence of fig. 5. The length of transition selective 27r pulse is identical as that 
of TT pulse but has double the pulse power, such that the dynamic phase remains refocused. Spectra corresponding 
to geometric phase shift gates (a) U[(3,4),(5,6)]=diag[l,l,-l,-l,-l,-l,l,l] implemented by applying (27r)y pulses on 
transitions C2 and C3. In this experiment the coherences hi, h2, and /ig acquire a Geometric phase of tt while the 
coherences ^3 and /14 acquire a phase of 27r. (b) U[(l,2),(3,4)]=diag[-l,-l,-l,-l,l,l,l,l] implemented by applying {2TT)y 
pulses on transitions Ci and C2 yielding the proton spectrum with coherences /13 and /14 yielding a tt phase shift, (c) 
U[(l, 2), (5, 6)]=diag[-l,-l, 1,1,-1, -1,1,1], implemented by applying (27r)y pulses on transitions Ci and C3, yielding the 
proton spectrum, in which all proton coherences acquire tt phase shift. 

(8) . Quantum circuit for implementation of two qubit Deutsch-Jozsa algorithm (Collins version). The pseudo 
Hadamard gate creates superposition of the four basis states. Uji is the unitary transformation, given in eqn.l4, 
corresponding to the function fi given in table (I). The final Hadamard gate in the original circuit is canceled by a 
detection pulse (psciido Hadamard gate) . Hence the result of the algorithm is available immediately after U fi . 

(9) . Schematic energy level diagram of dipolar coupled oriented CHsCN. The four levels of symmetric manifold 
are labeled as, basis states of a two qubit system. Transitions 2, 4 and 5 are degenerate, as shown in fig. 10(a). 

(10) . (a) Equilibrium proton spectrmn of oriented CH:>,CN. The labeling of transitions are in accordance with 
the notations used in fig. (9). The relative integrated intensities arc in the ratio 1.00:1.99:0.98 (theoretical ratios are 
1:2:1), and the splitting 3Dhh has a value 4968 Hz. (b) spectrum obtained by a small angle pulse after SOSS by the 
pulse sequence (tt)^ — — (tt)^ (fig. 11). The relative experimental integrated intensities arc in the ratio 3.02:3.98:2.97 
(theoretically expected 3:4:3). (c) spectnun obtained by a small angle pulse after the preparation of pseudo pure state 
|00), by using the pulse sequence, (tt)'^ — — Gz (fig. 11). 

(11) . The pulse sequence for the implementation of two qubit Deutsch-Jozsa algorithm (fig. 8), consists of four 
parts; (i) SOSS: (tt)^ — Gz — (tt)^, (u) Pseudo pure state |00): (tt)^ — (7r/2)^ — Gz, (in) Coherent superposition followed 
by f//7=diag[l -1 1 -1]: A {n/2)y on resonance pulse creates a superposition state = (l/2\/2)(|00) + \/3|01) + 
\/3|10) + 1 11)), Uf^ is implemented by [(27r)^ (27i')^^]. The duration of transition selective pulses (5.229 msec) is such 
that the evolution due to dipolar coupling makes a complete 2tt rotation and thus the phases are only due to cyclic 
geometric evolution of subspaces. (iv) Detection of single quantum coherences. 
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(12) . Implementation of two qubit Deutsch-Jozsa algorithm on oriented CH3CN, by using the pulse sequence 
shown in fig. 11. Transition selective pulses are of Gaussian shape with duration 5.229 msec. After pseudo Hadamard 
gate (fig. 11), controlled phase shift gate Uf. is implemented, (a) C//j is identity matrix which require no pulse, hence 
single quantum coherences are measured immediately after pseudo Hadamard gate. The integrated intensities are in 
the ratio, 1.69:3.00:1.72 (theoretically expected a/3:3:-\/3 ISl*!), confirms state {ipfi) of eqn.l5. (b) Uf^ (c) Uf^, are 
implemented by applying a single (27r)j, pulse on transitions 2 and 3 respectively, (d) Ufj is implemented by applying 
two (27r)j^ pulses on transitions 2 and 3. In fig. (a), none of the peaks are inverted, indicates the constant function. The 
fig.s(b,c and d) confirms the final states \1pf3), iV'/s) ^-nd \'ipf7) of eqn.(15), and atleast one of the peaks are inverted, 
indicates the balanced functions. 

(13) . Quantum circuit for implementing qubit-qutrit parity algorithm. Initial state is a mixed state followed by a 
pseudo Hadamard gate (h) on qubit register. Unitary operators of oracle are controlled phase shift gates, given in 
table(2) for various parities. The final step is the measurement on qubit register, which in NMR is the direct signal 
acquisition of resulting single quantum coherences. 

(14) . Schematic energy level diagram of oriented CH2FCN, consisting of symmetric and asymmetric manifolds. 
Flourinc transitions are labeled as Fi and F/ and that of protons as Hi. The transitions F2 and are degenerate 
(see fig. 15b). 

(15) . (a) and (b) are, respectively, equilibrium proton and fiourine spectra of oriented CH2FCN . In the proton 
spectrum, experimental integrated intensities are in the ratio 0.98:1.00:0.99:0.98 (theoretical ratios are 1:1:1:1), and 
that of fiourine are in the ratio 0.99:2.01:1.00 (theoretical ratios are 1:2:1). The transitions Hi and H^ (also H2 and 
Hi) are split by homonuclcar diploar coupling. The splitting is equal to 3Dhh which has a value 5694 Hz. The 
^^F-^H splitting in (b) is equal to {2Dfh+Jfh), which has a value 473 Hz. (c) Fiourine spectrum obtained by 
a small angle pulse after SOSS by the pulse sequence tt/2 — Gz ~ (Tr)^^'^'' — {Stt/A}^^ — Gz (fig. 16). The relative 
experimental integrated intensities are in the ratio 1.00:1.02:0.98 (theoretical ratios are 1:1:1). 

(16) . The pulse sequence to implement qubit-qutrit parity algorithm (fig. 13). The pulse sequence consists of four 
parts; (i) SOSS: n/2-Gz~{Tr)"^'"''-{3T:/4)'^^~Gz, (u) Pseudo Hadamard gate: On resonance {n/2)y'^ pulse, (in) Or- 
acle Of, shown in Table 2: During first t period three pairs of transition selective proton tt pulses [(tt)^^ i'^)fy+7r+7r/3)] 
[{'^)yH'^)fy+^+^/e)] [i^)y''i^)fy+7,-rr/2)] applied, which respectively implement cfj^ , (fj^ and d]j^J^ , and during 
second r period two pairs of transition selective fiourine tt pulses [('''')^^ ('''')(^^+7r+7r/6)] ('''')(^+7r+7r/3)] applied, 
which respectively implement c?J^^ and d^{^^ . The duration of transition selective tt pulses is set such that the evo- 
lution under the hetero nuclear couplings (J -I- 2D) makes a complete 2tt rotation, (iv) Detection of single quantum 
coherences of fiourine, from which the parity can be determined. 

(17) . Implementation of qubit-qutrit parity algorithm, by using the pulse sequence shown in fig. 16. The transition 
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selective tt pulses, for both proton and flourine, are of Gaussian shape with duration 12.85 msec. For odd (even) parity 
odd (even) number of peaks are inverted, in other words odd (even) number of single quantum coherences of flourine 
acquire a phase tt. (a), (b), (c) and (d) indicate the odd parity which are obtained by implementing Of, Of, Of and 
^3^'^'^^ (Table 2) respectively, (e), (f), (g) and (h) indicate the even parity, which are obtained by implementing I, 
0^^''^\ O^^'*^^ and O^^'^^ (Table 2) respectively. 
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